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Abstract 



CN| ' We present generalized and improved constructions for simulating quan- 

tum computers with a polynomial slowdown on lattices composed of qubits 
on which certain global versions of one- and two-qubit operations can be 
| performed. 

OO 

^ : 1 Introduction 

O 



In [9] we have shown that usual quantum circuits can be efficiently simulated 
on several regions of lattices composed of qubits on which instead of individual 
addressing of qubits we can perform certain global versions of one and two-qubit 
gates. The Hamiltonian of the global version of a two-qubit gate acting on a pair 
of a qubit is the (weighted) average of the translates of the Hamiltonian of the 
two-qubit gate. In [3] translations by lattice vectors were considered. A brief 
discussion and references to experimental results regarding implementation of 
global gates can also be found in . Our simulation requires a constant blowup 
of space and the slowdown is polynomial. 

Global two-qubit Hamiltonians considered in [5] and in this paper are spe- 
cial cases of Hamiltonians built from pairwise interactions. Results regarding 
efficient simulation of all such Hamiltonians using a fixed one with connected 
interaction graph and local one-qubit operations can be found e. g. in [5] (or 
in [T3] and [T2] in the context of NMR quantum computers). In [TT], the sets 
of operations which arc implcmcntablc as products of global gates from various 
collections are described. In the latter manuscript efficiency is not addressed: 
products of arbitrary length are allowed. 

Methods not requiring individual addressability of local spins are given in 
a series of papers by Benjamin and his collaborators. In [5] it is shown how 
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different leveis of addressability assumptions can be dispensed with in the case of 
a chain of spins. Multidimensional generalisations and issues related to physical 
implementations including error correction techniques are discussed in [TJ 31 [3] . 
In an intermediate proposal, given in [2J, spins only at even and odd positions 
have to be distinguished, but there is also a method that - like the proposals 
in [5] and the present paper - does not require that distinction either. Note 
however, that the endpoints of the chain play an essential role in the proposal 
given in [5J. In this paper symmetry receives a more attentive treatment and 
we do not make use of edge effects. The most closely related work to the 
present paper is probably [15] , where translation and reflection symmetric global 
operations in a chain of 5-level local systems were considered. 

In this paper we strengthen the results of [9] in two directions. First, we 
extend the method of [9] to more general models of global gates. In these models, 
the global Hamiltonians are built from translates of local ones by elements 
of groups more general than translation groups. In case of multidimensional 
lattices these groups may include not only translations but also transformations 
like rotations and reflections. Most notably, we propose constructions for the 
case of hypercubic lattices and global gates acting simultaneously on pairs of 
qubits at various Euclidean distances. In contrast to [S], the constructions in 
the present paper do not exploit that the domain containing the lattice points 
we are working with has borders. 

Another direction of improvement is restricting the set of global two-qubit 
gates needed for the simulation. We show that efficient simulation is possible 
with global one-qubit gates and global versions of two-qubit gates whose Hamil- 
tonians are diagonal in the computational basis. As the translates of such gates 
commute, in these cases the global gate is just the simultaneous action of the 
local translates. We note that some of the global gates in [3] are not of this 
form and therefore it seems to be much harder to implement them in practice. 

It is also natural to restrict the the " distances" of pairs of qubits which we 
let global gates act on. The most natural restriction would be allowing only 
next-neighbour interactions (e.g,, distance 1). Results with next-neighbor in- 
teractions can be found in [TJ] and [15]. The latter paper presents a method 
which works with next-neighbor interactions in lattices composed of 5-level sys- 
tems rather than qubits. The method of [14] works with a chain of qubits but 
it exploits that the chain has endpoints. It is not known if there is a version 
of our method using only global two-qubit gates with distance 1. In the one- 
dimensional case we can present a version in which the global two-qubit gates 
have distance at most 22. It works with qubits and does not make use of the 
border of the chain. Furthermore, our method also works if there is certain 
imperfection in the global gates. 

The constructions of the present paper require an initial state where qubits 
at certain positions are set to one while the others are set to zero. Obviously, 
form the all-zero state such state cannot be achieved using global gates. In [15] 
this problem is circumvented using certain techniques such as running several 
reflection symmetric simulations at sufficiently large distances simultaneously. 
Such techniques seem to be applicable in the case of 2-level (qubit) systems as 
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well - maybe using a bigger system of global two-qubit gates in the initialization 
phase. 

The structure of this paper is the following. The model of global gates 
is described in Section [21 Our main theorem is stated in Section [3J where 
also combinatorial notions describing the schemes making efficient simulation 
possible are introduced. Section \E\ is devoted to examples for such schemes, 
while the proof of our main result can be found in Section [4j We present a 
construction for simulation using only global two-qubit gates with distances 
between 1 and 22 in Section [5) 

2 Global gates 

In this section we introduce some notation and describe a general mathematical 
model for global one- and two-qubit gates. 

Throughout the paper we assume that G is a transitive permutation group 
on a possibly infinite set and D is a subset of ft. For p <E 51 and g € G we 
denote by p 9 the image of p under g. The permutation action of G induces an 
equivalence relation ~ on D x D: (p, q) ~ (p' , q') if there is an element g € G 
such that p' = p 9 and q' = q 9 . We denote by O the set of the equivalence 
classes of ~ different from the diagonal {(p,p)\p & D}. For (p, (?) € D X D the 
equivalence class containing (p, q) is denoted by C( p q y 

We encourage the reader to consider the following instructive class of exam- 
ples. Here Vt = R s , the s-dimensional Euclidean space and G is the Euclidean 
group E(s) which consists of isometries of R s . For every subset D C M. s , two 
pairs (p, q), (p', q') ofDxD are equivalent under the action of E{s) if and only 
if the distances \p — q\ and \p' — q'\ are the same. As a consequence, the classes 
of O can be indexed by positive real numbers (distances) 8 and (p, q) £ Cg if and 
only if \p — q\ = S. (Of course, depending on the choice of D, the class Cs may 
be empty for certain 5's. We shall refer to this geometric example in intuitive 
explanations of abstract notions and arguments. 

To define global gates, we assume that D is finite with \D\ = n and we have 
a configuration of n qubits sitting at each element of D. The Hilbert space of 
the pure states over these n qubits is C 2 . The elements of the standard basis 
are indexed by the functions a : D — > {0, 1}. In order to shorten notation, for 
p 6 D we also write a p for the value a(p) e {0, 1}. 

As a part of our model we introduce a balance function W : D x D — * K >0 
on the pairs of D. This function will make it possible to model effects like 
imperfection of global gates, namely that the tool performing a global gate (e.g, 
a laser beam or a magnetic field) may act with different strength on pairs or 
singletons of qubits at different positions. Of course we could take W to be 
constant. This choice would express invariance of global gates under the action 
of G perfectly. In order to shorten notation we will also denote W(p,p) by 
W(p). 

Below we give a formal definition of the model of global gates. The first 
part is devoted to describe how global Hamiltonians are built from local ones. 
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Then global gates are obtained from global Hamiltonians in the usual way. 
As in Section [4] we shall make use of the formalism introduced here in the 
context of arbitrary operations, in the first part we do not assume Hermiticity of 
operations. That is, by an operation we mean an arbitrary linear transformation 
or matrix of the appropriate dimension. However, the reader not interested in 
the details given in Section [?] may think merely of Hamiltonians first. 

For a one-qubit operation or 2 x 2 matrix M, whose rows and columns are 
indexed by and 1, and an element p G D we write M p for the rt-qubit operation 
which acts as M on the qubit at position p: 

M p = i Ma "- b p if ° s = bs for every s e D Mp}> 
a ' b [0 otherwise. 

We will refer to M p as a local one-qubit operation at position p. The global 
one-qubit operation corresponding to M is 

M° = W(p)M p . 

p£D 

Similarly, for a 2-qubit operation or 2 2 x 2 2 matrix M (rows and columns indexed 
by 00, 01, 10, 11) and a pair of elements p ^= q G D, the local two-qubit operation 
M( p ' q } at position (p, q) is defined as 

m (p,9) = f M (a p ,a q ),(b p ,b q ) if a s = b s for every s G D \ {p, q}, 
"■■ b \ otherwise, 

and for a C G O the global two-qubit operation is 

M c = J2 W(p,q)M ( P' q l 
(p,«)ec 

A two-qubit global Hamiltonian is a matrix of the form H c where H is an 
Hcrmitian 2-qubit operation and a global 2-qubit gate is an operation of the 
form exp(— iH c ) where H c is a global two-qubit Hamiltonian. Global one-qubit 
gates are of the form exp(— iH°) where H is an Hermitian one-qubit operation. 

Note that global one-qubit gates can be interpreted as parallel execution of 
local operations since 

exp{-iH°) = Yl cxp(-iW(p)H p ), 

where the product on the right hand side can be taken in an arbitrary order 
of the terms. The analogous statement for global two-qubit operations is not 
true in general. Only if for all (p, q), (p' , q') G C the local operators H^ p,q ^ and 
ff(p ,i ) commute (this happens for example if H is diagonal in the standard 
basis) can exp(— iH c ) be decomposed as Y[, p q ^ eC exp(—iW(p, q)H^ P:q ^). 
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3 Simulation results 



An intuitive description of our method for simulating a " usual" quantum com- 
putation with global gates acting on D is the following. We will designate two 
disjoint subsets P and R of D. The qubits at positions in R will serve as refer- 
ence points. Before and after each logical step of the simulation, their value will 
be set to one. Similarly, the values of the qubits at positions outside P and R 
will be set to zero. The subset P will be the workspace. The qubits at positions 
in P can take arbitrary values and in each logical step we simulate a local one- 
or two-qubit gate at certain (pairs of) positions in P. To be more precise, each 
logical step will approximate such a gate. 

By approximation we mean approximation in terms of the operator norm. 
For an operator U on the Hilbert space C 2 we denote by \\U\\ the operator 
norm of U: \\U\\ = sup^i^ \Ux\. Note that ||AB|| < \\A\\ ■ \\B\\. Furthermore, if 
||j4i||, . . . , ||Ajv||, ll-Bill, . . . , \\Bn\\ < 1 (this holds in particular if Ai and B>i are 
unitary operators), then we have 

N 

\\Ai---A N -Bi---B N \\ <^2\\Aj -Bj\\. (1) 

i=i 

We say that A e-approximates B if \\A — B\\ < e. Equation |T]) implies that 
in order to e- approximate a circuit consisting of t gates it is sufficient to e/£- 
approximate each gate occurring in the circuit. 

Below is a very informal - metaphoric - description of our simulation method. 
In this description qubits correspond to two-state ("up" or "down") switches 
placed at positions in D. The global one-qubit gates are be modelled by the 
ability of flipping the state of the switches simultaneously. The additional device 
corresponding to the two-qubit gates with distance S is the following. For every 
p £ D, this tool commits the change of switch at p if there is a position q at 
distance S such that the original state of the switch at q was " up" . The value of 
the commitment is defined as the number of such points q. After a sequence of 
commitment operations the change is performed at positions where it has been 
committed by all of the commitments. 

We work with states where the switches at points in R will be always " up" , 
switches at points in P can be both " up" and " down" , while the rest are always 
" down" . The goal is to be able switching the state of a specific switch in p e P 
by keeping the state of the others using the operations described above. We also 
assume that this must be done performing the global flip operation followed by 
a sequence of commitments which depend only on the position p and not on 
the particular state of the switches in P. Furthermore, we require that for 
committed flips, the total value of commitments does not depend either on the 
particular state. To perform this we need switches at r%, . . . , in R such that 
in any case, commitment operations with distances \p — r*i|, . . . , \p — commit 
the flip only at position p. That is, we need to be able to uniquely locate point 
p e P in such a strong sense just by the distances of it form points of R. We 
shall use the term "addressable" to express this property. 
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We stress that the description above is only metaphoric, looking for any 
direct connection with the reality would be misleading. The actual and accu- 
rate description of the simulation method can be found in the rather technical 
Section 0] Below we give the formal definitions. 

We call a pair (P, P) of disjoint subsets of D a workspace-base scheme. We 
say that the workspace-base scheme (P, R) admits the function a : D — > {0, 1} 
(or a is admissible if (P, R) is clear from the context) if a r = 1 for every r 6 R 
and a q — for every q e D \ (P U R). We say that p € P is addressable 



by T-i, . 



, a- 



G P if for every fc + 1-tuple (p' , 



y k ) e D x (PU P) fc , the 



property that for each i = 1, . . . , k either (j>' , r'A £ C( p , ri ) or {r J i ,j/) £ C{p,n) f° r 
each i — 1, . . . , k implies p' — p and r[ . . . , r' k € P. We say that p e P is fc- 
addressable if there exist ri, . . . ,r k € P such that p is addressable by ri, . . . , r^. 
We call (P, P) fc-addressable if every p € P is fc-addressable. 

The subspace of C 2 spanned by the basis elements corresponding to the 
admissible functions is obviously isomorphic to C 2 ' P ' . Its elements (of length 1) 
are called admissible vectors (or states, respectively). The remaining basis ele- 
ments will span the orthogonal complement of the subspace of admissible states. 
The vectors (or states) in the latter subspace will be called inadmissible. Our 
first (and main) goal is to produce unitary operations which (approximately) 
preserve the subspace of admissible vectors and, restricted to this subspace, (ap- 
proximately) act as local one-qubit gates. In other words, if we adopt an order 
of the basis where the first 2' p l basis elements correspond to the admissible 
functions and the rest correspond to the inadmissible functions, we will produce 
operations whose matrices are (approximately) block diagonal and have upper 
left 2' p l x 2' p l block which is (approximately) the same as the matrix of a local 
one-qubit gate acting on a qubit at position peP. Then, having the local one- 
qubit gates at hand in order to get a universal quantum computer on on the 
admissible states, it is sufficient to use a unitary operator which has an all-pair 
entangling Hamiltonian on the admissible states. 

We shall achieve our main goal by taking commutators of certain global one- 
qubit Hamiltonians with Hamiltonians of the form A c at the Lie algebra level, 
where C G O and A is the Hamiltonian of the so-called controlled phase shift 
operation. The entries of the matrix of this two-qubit Hamiltonian are all zero 
except in the lower right corner (corresponding to the basis element indexed by 
11): 



A = 



The corresponding global two-qubit gates are 



( 
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(2) 



exp 



/ / 





\ \ 







-T, J 



\ C (p,r)\ 



, (T eM.p € P,r € P), 



(3) 



G 



where i stands for the imaginary unit \f— T. 

We shall approximate the local one-qubit gates by products of global two-qubit 
gates of the form and global one-qubit gates of the form 



exp 







(z e C). 



(4) 



In Section 2] we shall prove the following. 



W(p',q') 



< W. 



Theorem 1 Assume that for every p',p",(/,q" G D we have w ^ p „ q ,,\ 
Assume further that (P,R) is a workspace-base scheme on D. Then, for every 
k-addressable p € P, and for every < e < 1, every operation which acts on 
admissible states as a one-qubit gate at p with Hamiltonian of norm bounded 
by a constant can be e- approximated by a product of (nw/e)°( k > global gates of 
type 0j and Q). 

If the "distance" between two points p,q S P cannot occur as a "distance" 
between a point p' in P and another point r in R, that is 

(Cp, g U C q}P ) n C p '^ r — for every p, q,p' € P and r £ R, (5) 

then the global gate 



exp 



/ / 





\ \ -Ti ) 



(6) 



/ 



acts on the admissible states as a unitary operator whose Hamiltonian entan- 
gles the qubits at p and q (and also the pairs of qubits in P with the same 
"distance"). Global gates of this type, together with the simulated one-qubit 
gates on admissible states, using the selective decoupling techniques (see e.g. 
p~2\ IT3l 15] ) . give an efficient simulation of any quantum circuit. We obtain the 
following. I 1 ] 

Corollary 2 Assume that for every for every p' , p" , q' , q" G D we have w(p" 'q") — 
w. Assume further that (P, R) is a k-addressable workspace-base scheme on D 
with additional property (0). Then, for every < e < 1, the action of any 
quantum circuit of length £ on P on admissible states can be e- approximated by 
a product of (nw£/e)°( k ) global gates of type (0), Q and (0. 

Here by a quantum circuit we mean a circuit built from 2-qubit gates, that 
is a sequence of unitary 2-qubit operations. As 2-qubit gates are universal (see 
e-g- [Z])) any circuit in the sense of can be efficiently simulated as well. 



1 We are indebted to an anonymous referee for suggesting the argument above. This has 
made it possible to replace a much more complicated argument - and a somewhat weaker 
statement - in a preliminary version of this paper. 
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4 Extracting local gates by taking commutators 



Our main tool will be taking commutators with global two-qubit Hamiltionians 
of the form A c , where the matrix A is as in Q . This will be the tool analogous 
to the commitment operation in the metaphoric description. See Proposition [7] 
for the effect of an appropriate sequence of such commutations. It is easy to see 
that, for every C G O, A c is a diagonal matrix whose entry at position a, a is 



A„ 



E W M- 



(p. St) £ C 



It follows that for arbitrary a,b : D — > {0, 1}, 



( 



\ 



(p. g) e c 

bp = b a = 1 



where i? aj b is the 2> D > x 2'"°' elementary matrix having one at position a, b and 
zero elsewhere. For a : D — > {0, 1} we denote by d p a the function which takes 
value zero at p and coincides with a on D \ {p}. That is, d p a is obtained from 

a by deleting the p^ bit. Commutators of A c with elementary matrices of the 
form E ay Q pa turn out to be of particular interest. We have [A c , E a ^g pa ] = if 
a p = and 



[A c ,E a , dpa ] = 



£ W(p,q)+ £ WfoP) 



<p> <j) e c 

a„ = 1 



(<!, P) 6 C 
a„ = 1 



^ Q ,a DQ (7) 



if a r 



by 



For C € C we define a modified balance function Wc on Dx D\{(p 1 p)\p G D} 



W C (p,q) 



if (p,q) and (q,p) g" C, 

W(p, q) if (p, ?)eC and (g,p) g C, 

if (g,p) eC and (p,g) 0C, 
W(p, q) + W(q, p) if (p, q) eC and (g, p) G C. 



Iterated application of ([7|) gives the following two lemmas. 

Lemma 3 Let p G D, Ci, . . . ,Ck G O, and a : D —> {0, 1} with a p — 1. Assume 
further that there is no k-tuple q^, . . . , q' k G D such that a q ^ — . . . = 
either (p, q[) G C; or (g£,p) G Ci holds for every i — 1, . . . . fc. 27ien 

[A c M...,[^,S a ,a pa ]]] = 



a„t = 1 and 
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Lemma 4 Let p € D, C\, . . . ,C& € O, and a : D — ► {0, 1} wz£/i a p = 1. Let 
Q C Z? \ {p} smc/i i/iaf a 9 = 1 /or euen/ g G Q and qi, . . . ,q k G Q swc/i i/ia£ 
(P) 9i) € Ci or {<liiP) G Ci / or i = 1) • • • j k. Assume further that if for a k-tuple 
q[, ■ ■ ■ , q' k G D we have a g j = . . . = = 1 and either (p, q^) S C, or (g^,p) G Ci 
holds for every i = 1, . . . , k, then q[, . . . , q' k G Q. Then 

[A Cl ,[..., [A c \ E a<9pa ]]] =Wx--- W k E a , dpa , 

where W l = J2 q eQ W d (P) «)■ 

Lemmas S] and [3] imply the following. 

Lemma 5 Let (P, P) oe £/ie workspace-base scheme, a : D — > {0, 1}, p G P 
and q €z D. Assume that p is addressable by r%,...,rk G P. Under these 
assumptions, if a is admissible then 

\ A Ci r mc* P in- / " ' ' w hE a ,d q a if 'q=P and Op = 1, 
[A ,[..., [A , E a , dqa \\\ - | o otherwise, 

where Ci = Ct p , ri ) an d Wi = X^rei? i^c. (p, r ) f or i = l,..-,k. If a is not 
admissible but d q a is admissible then 

[A c \[...,[A c ",E a , dqa ]]}=0. 

Proof. Assume that either a or d q a is admissible and that the commutator is 
nonzero. Then a q = 1, and, by Lemma there exist r[ . . . ,r' k G D such that 
a.,.j = . . . = a, r ^ = 1 and either {q,r[) G Ci or (r^,g) G Cj. Since a or 9 9 a is 
admissible, r\ G PUR and from the definition of workspace-base schemes it 
follows that q—p and r[ = Ti for i = 1, . . . , fe. This shows the second statement 
and part of the first one. The rest follows from Lemma 2] □ 

For an elementary one-qubit global operation we have the following. 

Lemma 6 Let (P, R) be a workspace-base scheme, a, 6 : D — > {0,1}, and as- 
sume that p € P is addressable by r\, . . . ,r k . Suppose further that either a or b 
is admissible. Then, for the 2x2 matrix 



B 



1 




,6' 



we have 

[A c \[...,[A c \Bl b ]]]=W(p)W 1 ---W k Bl 
where d = C( p , ri ) and Wi = J2 r eR W ^ (Pi r) for i = 1, . . . ,k. 

Proof. We denote by C the commutator on the left hand side of the asserted 
equality. For every ijeDwe have B q — J2 a \a q =i E a ,d q a- Therefore 

b* = J2 W ^) E *w=E E w (i) E ^ qa 

q&D a\a q = l a q \a q = l 
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From this equality, using Lemma O we infer 

c = f W(p) IlLi W (P> r i) if b = and fl p = !' 

a ' \ otherwise, 

whenever either a or & is admissible. From the latter equality the assertion 
follows because B v a b is one if a v = 1, b = d p a; and zero otherwise. □ 

For a global one-qubit Hamiltonian we obtain the following. 

Proposition 7 Assume that (P, R) is a workspace-base scheme on D, and let 



X 



z 
-z 



where z is a complex number. Suppose further that p £ P is addressable by 
ri, . . . , rfc G P. TTien /or every pair a,b : D —> {0, 1} suc/i that either a or b is 
admissible, we have 

HA Cl ,[..., [-iA c - , X*]]] a , h = WCpjW! • • • w k xi b , 

where Ci = C (?vr .) 7 = ^r-Gfl ^ (j>, r), and 



X = 



{-i) k z 
-? fc z 



Proof. Observe that X = zB — zB> where B is the same 2x2 matrix as in 
Lemma [6] Hence, using also that the matrices A Ci are self-adjoint, we obtain 

[-iA c i ,[..., {-tA c » , X-}}} = (-i) k [A e > ,[..., [A c * , X-}}} = 
z(-t) k L4 C > ,[..., L4 C * , B*]]] - [A Cl ,[..., [A c * , B*t]]] = 

z(-z) fe L4 Cl ,[..., [A c * , -B°]]] — z(~z) k [A c - ,[..., [A c * , fl<>]]]t 

Using Lemma [31 this equality gives 

HA c M...,HA^,^]]] Q , b = zH^w^-^b^ 

-z(-i) k W{p)W 1 ---W k Bl P b 

whenever either a or & is admissible. □ 

Proposition [7] is used to prove Theorem [TJ i.e., to show that the local one- 
qubit gates can be efficiently approximated. We need the following fact on 
approximation of an operator whose Hamiltonian is a commutator. For the 
proof, see e.g., [S]. 
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Fact 8 There is an absolute constant c > 0, such that 

N 

exp(--^X-) • exp(--^y-) • exp(--^X) . exp ( --^Y) 



exp([— iX, —iY] 



< c ■ M 3 N~? 



for any N > M 2 , where X and Y are Hermitian operators on the Hilbert space 
C 2 " and M = max{||X||, ||Y||,1}. 

We remark that commutators (both group and Lie theoretic, in view of 
Fact[8j also play an important role in proof the Kitaev-Solovay theorem (see e.g., 
|10j). It turns out that decomposing unitaries as group theoretic commutators 
can be used to exponentially improving a sufficiently good approximation of a 
unitary operator by a circuit built from a given gate set. 

Now we are in a position to prove Theorem [1] 

Proof of Theorem Q3 We show the statement for gates of type 

exp(Y _°_ I jy (zeC,|z|<l). (8) 

Then, the general case for a general one-qubit gate follows from Fact [5] and 
Trotter's formula since the third basis element 

' i 
-i 

of the Lie algebra SU2 equals ^ of the commutator of the first two ones which 
are 

ur v = ( _° : \ ) and ^ K = ( ° I 

Let 

r, ( z 



W(p) \ -z 

and for i = 1, . . . , k let X t = l/W t A c [p r%) where W t = EreR ^V^, (p, r). Then, 
by Proposition [7J on admissible states the operator exp([iXi, [. . . [— iXk, —iZ]]]) 
acts in the same way as the required one-qubit local gate of type ((5]). For 
i = 1, ...,k we denote the commutator [iXi, [. . . [— iXk, — iZ]]]) by Yi. Also 
set Yfc + i = — iZ . Then for i = 1, . . . , k we have Yi — [— iXi, li+i]. We have 
\\Z\\ = 0(nw), \\Xi\\ = O(nw) and = 0({2nw) k - z+1 ). It follows that all the 
norms (i = 2, . . . , k+1) \\Xi\\ (i = l,...,k) are bounded by M = 0(2nw) k . 
Set ei = e. By Fact [SI we obtain an ei/2-approximation of expY^ by a product 
containing 2Ni terms of exp(±Xi) and 2iVj terms of exp(±Yi + i), where Ni — 
0(M 6 )/ei. If we substitute each factor exp(±Vi + i) by approximations with 
error ej+i = Ci/4iVj, the total error of the product will be bounded by ej. For 
the sequence Q we obtain a recursion e^+i = 0(M 6 £i) (i = 1, . . . , k — 1). That 
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is, a = (c'M 6 ) 1 -^ for some constant d and hence Ni = 0((c"Af 6 )'). The total 
number of terms is 47Vi + 4iV x ■ AN 2 + . . .+47Vx • • • AN k = 0(k(4c" M 6 ) k( - k+1 ^ 2 ) = 



5 Addressable workspace— base schemes 

In this section we give some examples of addressable workspace-base schemes 
corresponding to various group actions. In order to obtain nice ("uniform") 
constructions, we extend the notion of addressability (and connectedness) to 
possibly infinite subsets D C VL in the obvious way. 

Assume that (P, R) is a workspace-base scheme in D with R = {ri, . . . , r^}. 
We say that (P, R) is strictly addressable if for every p £ P and for every k + 1- 
tuple (p',r[, ...,r' k )£ Dx(PUR) k (p',r-) £ C (P:n )UC (nj) ) (i = l,...,k) implies 
p' = p and r\ = r j (i = 1, . . . , k). This means that every p £ P is addressable by 
the whole of P in a strict sense where - in terms of our geometric example - the 
sequence of "distances" of p from r% identifies p and r±,...,rk- We note that 
strict addressability implies condition §5§ and therefore Corollary [2] is applicable 
to strictly addressable schemes. 

5.1 Schemes for translation-invariant global gates 

This part is devoted to constructions of strictly addressable workspace-base 
schemes for translation- invariant global gates in abelian groups, especially in 
the s-dimensional lattice Z s . In particular, we shall show that (Pi,R±) with 
Ri = {1,3,9} and Pi = 4Z = {...,-12,-8,-4,0,4,8,12,...} is a strictly 
addressable scheme for translations in Z of density 1/4. 

The model of translations in Z can be interpreted as that the qubits sit in 
a doubly infinite chain and the group acting on the chain is given by the shifts. 
Assuming a constant balance function, the global gates are shift-invariant. Just 
like in Section^ the actual (finite) version of the model is obtained by restricting 
the action of the global operation to a finite section of the chain. 

For all positive integers s the scheme above lifts to strictly addressable 
schemes (P S ,P S ) where R s = {(1, 0, . . . , 0), (3, 0, . . . , 0), (9, . . . , 0)} and P s = 
{(4zi, Z2, ■ ■ ■ , z s )\zi £ Z}. It is an open question whether 1/4 is indeed the opti- 
mal density or there are more efficient constructions. We start with the formal 
definitions and give some lemmas which may be useful in constructing other 
addressable schemes. 

Let G be an abelian group. We use the additive notation for the group 
operations. The group G acts on ft = G by shifting: p 9 = p + g. We set D = G. 
As (p, q) ~ (p 1 , q') iff p' — q' = p — q, the elements of O can be indexed by the 
elements of G \ {0}: (p, q) £ C g iff p — q = g. In other words, Ci p , q ) — C p - q . 

Under some mild restrictions, we can lift strictly addressable schemes from 
factor groups as follows. 

Lemma 9 Assume that <j> is a surjective homomorphism from the abelian group 
G onto K , and (Q, S — {si, . . . , Sk}) is a strictly addressable scheme on K . Let 
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P = <f)~ 1 Q. For each i = 1, . . . , k we pick a single element ri from {sj}. 
Then, assuming that there exist i,j G {1, . . . k} such that 2r^ ^ 2rj, the scheme 
(P, R = {fx, . . • , rfc}) is a strictly addressable on G. 

Proof. Assume that p' - r[ = ±(p - n) with p G P, p' G G, r< G P U R (i = 
1, . . . , k). Then 0(p')-^(r-) = ±(<p(p) - <p(n)) and ^(p) G Q and 0(r{) G QUS, 
and strict addressability of (Q,S) implies 4>(p') = 4>{p) and 4>( r i) — <?K r f° r 
i = 1, . . . , k. From the latter facts and disjointness of Q and S we infer that 
r\ = T{ for i = 1, . . . , k. Hence p 1 ' — = ±(p — rj) for i = 1, . . . , k. Yip ^ p' then 
this is only possible if 2r^ = p + p' for every i. □ 

Lemma Ogives several straightforward ways of constructing strictly address- 
able schemes in multidimensional lattices from one-dimensional schemes. In 
particular, strict addressability of the scheme (Pi, Pi) given above implies that 
of (P s , R s ) for every s. 

As another application of Lemma [9] we can give a simple construction for a 
strictly addressable scheme of density 1/14 in Z. 

Proposition 10 In Zi4, the additive group of integers modulo 14, the scheme 
{0},{1,3,7} is a strictly addressable workspace-base scheme. As a consequence, 
so is the scheme {\Az\z € Z},{1,3,7} in Z. 

Proof. Assume that for p' G Z14 and for r[ , r' 2 , r' 3 G {0, 1, 3, 7} we havep' — r[ = 
±1, p' — r' 2 = ±3, p' — ?'3 = ±7. If p' is odd then r[ , r' 2 , r' 3 must be even therefore 
r[ = r' 2 — 73 = and hence {±(p' — r[)\i = 1,2,3} = {±p}, a set containing 
at most 2j3 elements, a contradiction. Therefore p' is even and r'^r'^r'^ must 
be odd which is only possible if {r[, r' 2 , r' 3 } C {1,3,7}. On the other hand, 
{±(2 - 1), ±(2 - 3), ±(2 - 7)} = {1, 13, 9, 5}, {±(4 - 1), ±(4 - 3), ±(4 - 7)} = 
{3, 11, 1,13}, {±(6-1), ±(6-3), ±(6-7)} = {5, 9, 3, 11, 13,1}, {±(8-1), ±(8- 
3), ±(8 -7)} = {7,5,9,1,13}, {±(10 - 1), ±(10 - 3), ±(10 - 7)} = {9,5,7,3} 
and none of these sets contain {1, 3, 7}, therefore p' — 0. As —1 and —3 are not 
in {1,3,7}, the only possibility is that r\ = 1, r% = 3 and r 3 — 7. The second 
statement follows from Lemma |9] □ 

We proceed with the somewhat tedious proof of strict addressability of the 
scheme (Pi, Pi) defined above. In order to support similar but potentially more 
efficient constructions, we give the proof in form of some statements. 

Note that strictly addressable schemes (P, R) where P is non-empty and 
|P| < 2 do not exist in any abelian group G. In the remainder of this subsection 
we consider schemes with three-element bases. 

Lemma 11 Let G be an abelian group, ri,r 2 ,r 3 G G such that r^ — rj ^ Ty—Ty 
for {i,j) 7^ (i', ]')■ Then for every p G G such that 2p ^ {ri + r 2l 7*1 + r 3 , r 2 + r$} 
for every p' G G, for every r' x ,r' 2 , r' 3 G {r±, r 2l r 3 } the condition p—ri = ±(p'—r / i ) 
(i = 1, 2, 3J implies p' = p, r[ = T\, r' 2 — r 2 and r' 3 = r 3 . 
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Proof. As 2p g" {ri + r 2 , r% + r 3 , r 2 + r 3 }, we have p — ri =/= ±(p — rj) for i ^= j. 
Therefore {r[, r 2 , r' 3 } — {ri, r 2 , r 3 } and there exists a permutation tt of {1, 2, 3} 
and v\, i/2, v 3 S {±1} such that p — ri = i/i(p' — r^). There exist i ^ j such that 
Vi = i/j = v. We have p — Ti — v{p' — r,-« ) and p — rj = v{p' — rj*). Subtracting 
these equalities we obtain Tj — r% — v{rj* — rj») and hence {i,j} — 
As a consequence for the third element k S {1,2,3} \ {i,j} we have k 7 ' = k. If 
p' =/= p then p' — rfc = rfc — p (that is, p' = Irk — p) and regarding p' — ri and 
p — rj we are left with the following cases. 

(1) v = — 1 and 7r = 1 is the identity 

(2) v = 1 and tt = (12) 

(3) 1/ = -1 and vr = (12). 

In case (3) we have p— rj = r j — p' and p — rj = n—j/, that is, p' = r,; + rj — p. 
Comparing this with p' = 2r& — p we obtain 2rfe = r^ + rj and hence r^ — r^ = 
rfc — rj, a contradiction with our assumptions. 

In case (1) we obtain p' = 2r^ — p = 2rj — p = Irk — P, which is solvable for 
every p iff 2r^ = 2rj = 2r^. 

In case (2) we obtain p' = p — rj + rj = p — rj + r^ = p — 2r& which is solvable 
iff 2ri = 2rj and 2p — 2r^ — ri + rj. □ 

Proposition 12 Let G be an abelian group and assume that there is an epi- 
morphism tj> : G — > Z 4 . We set D = G, R = {ri,r 2 ,r 3 } and P = (/> _1 (0) 
G\2p — ri+rj for some i,j = 1,2,3}, where <fi(R) = {1,3}, andr.i — rj ^ ry — rji 
whenever i ^ j and ^ (i',j')- Then (P, R) is a strictly addressable scheme 
in G. 

Proof. Let p E P, p' € G, r' l7 r 2 ,r 3 e P U R such that p' - r[ = ±(p - n). 
Observe that {<t>(p-n)\i = 1,2,3} = {1,3}. Therefore {cf>(p' - n)\i = 1,2,3} = 
{1, 3} as well. As 4>({P U {n, r 2 , r 3 }) = {0, 1, 3}, this is only possible if 4>{p') is 
even and </>(r-) are odd (i — 1,2,3). As P C _1 (O), this implies {ri,r 2 ,r 3 } C 
{ri, r 2 , r 3 } whence, by Lemma [TTl p' = p and r^ = r^ (i = 1, 2,3). □ 

Proposition [T2l gives that (Pi, Pi) is a strictly addressable scheme in the 
group Z, where i?i = {1,3,9} and Pi = 4Z. In Z 4 m for to > 2 we obtain the 
strictly addressable scheme (Z4 m \ {2m + 2, 2to + 6}, {1, 3, 9}). 

We shall make use of the following lemma in Section [6] 

Lemma 13 Let G be an abelian group and assume that there is an epimor- 
phism cf> : G — ► Z4. Let D = G and let ri,r 2 ,r 3 ,r4 £ (f)^ 1 ^,^} such that 
4>({ri, r 2 , r 3 }) = {1,3}, 3r 4 7^ ri + r 2 + r 3 , and r^ — rj 7^ r^ — rj> whenever 
i ^ j and 7^ (i',f). Then with P = _1 (O), the element r 4 is strictly 

addressable by {ri,r 2 ,r 3 } with respect to the scheme (P U {r 4 }, {ri, r 2 , r 3 }). 

Proof. Assume that for some r 4 6 G, for r^r^rg 6 PU {ri, r 2 , r 3 , r 4 } we 
have r 4 - r^ = ±(r 4 - r,) (i = 1,2,3). We have {0(r 4 - r 4 )|i = 1,2,3} = 
{0,2}, whence also {0(r 4 - r' t )\i = 1,2,3} = {0,2}. Since <jf>(P) = {0} and 
<t>{ r ii r 2, r3, r 4 } = {1, 3} this implies that </>(r 4 ) is odd and r[, r 2 , r 3 g {ri, r 2 , r 3 }. 
There exist /1 = ±1 and i =/= j G {1, 2, 3} such that r' A — r[ = /i(r 4 — ri) and 
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r' 4 — Tj — /i(r4 — rj). As a consequence, r[ — r'j = ju(rj — rj), whence either 
jit = 1, = Ti, Tj = rj or /i = — 1, = rj, r^ = rj. In the first case we have 
r' 4 = r4 and, by the properties of {ri, r2, r"3, r^}, also rj. = rj, for fc = 1,2,3. 
Assume the second case. Then r' 4 = r^ + rj — r4. Let k = {1, 2, 3} \ {i, j}. Then 
ri + rj —T^ — r' k — r' A — r' k = v{r±~Tk) for some v = ±1. If 1/ = — 1 then we obtain 
r^ — rj, = rfe — rj, a contradiction to the properties of {n, r2, r3, r4}. If 1/ = 1 
then 2r4 + rj-. = r^ + rj + r^. If r' k = r^, then we obtain r4 — r^ = — ?'4, while 
if rl = r4 then 3r4 = rj + rj + rfc = ri + r2 + r3. Both possibilities contradict 
to the properties of {ri, r2, r3, r4}. □ 

5.2 Schemes for Euclidean global gates in If 

In this subsection we show that strictly addressable schemes (P, R) exist in the 
lattice Z s where \R\ = min{3, s + 1} and P has density 1/8 S in Z s . Here 51 = K s , 
D = Z s , G = E(s) and the orbits of G on pairs correspond to distances. This is 
of the type of our main geometric example for global gates given in Section [2J 
For s — 1 see Proposition [T2l gives a scheme of density 4. The optimality of the 
density (l/8) s (for s > 1) is an open question. Also note that the addressable 
schemes for the group E(s) are also addressable schemes for the subgroups of 
E(s), most notably for the the group of isometries preserving the lattice Z s . 

Assume that s > 1. We shall actually show that for an appropriate choice of 
an s + 1-clement subset R of Z s , and a set P consisting of almost all points of 
the set 8Z S = {(8zi, . . . , 8z s )\zi € Z}, the scheme (P, R) is strictly addressable. 
Let k — s + 2 and for i = 1, . . . , k write Pi for the vector of variables Pi — 

{xi\ , Xi2 , . . . , Xi s ) . 

We set k = s + 2. A /c-element subset S C K s is reconstructible, if any other 
/c-element subset S" C K s which has the same distribution of pairwise distances 
among its points is actually isomorphic to S with respect to rigid motions (i.e., 
the group E(s)). 

From a result of Boutin and Kemper (Theorem 2.6 in [5]) it follows that 
there exists a nonzero polynomial /(Pi, . . . , Pfc) with real coefficients such that 
if pi,p2, ■ ■ ■ ,Pk € R s are points with /(pi, . • • ,Pfe) 7^ then {pi,...,pk} is 
reconstructible. The polynomial / has the additional property (Lemma 2.9, 
loc. cit.) that /(pi, . . . ,pk) = if the pairwise distances among the pi are not 
all distinct. 

If ^ g(x) €E WL[x] is a polynomial and H is an infinite subset of R, then 
there exists a c € P such that g(c) ^ 0. A repeated application of this simple 
fact to the Boutin-Kemper polynomial shows the existence of a subset R = 
{ri, . . . , r^^i} C 4Z S + (1, 0, . . . , 0) such that the polynomial 

h(P k ) :=f(ri,r 2 ,...,r k - 1 ,P k ) 

is not identically zero. Moreover, we can suppose that v\ € 8Z S + (1, 0, . . . , 0) 
and r2 £ 8Z S + (5, 0, . . . , 0). Let P be the set of points p from 8Z S for which 
M.p) 7^ 0- Note that P contains almost every point of 8Z S . Also, if p 6 P 
then the s + 1 distances \p — 7*4] are all different, and the set of these distances 
uniquely determines p among the points of R s . 
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We now show that the above pair (R, P) is a strictly addressable scheme on 
D = Z s . We note first that for every p £ P we have \p— ri\ 2 = 1 (mod 16) while 
\p — t2\ 2 = 9 (mod 16). Assume that we have a p £ P. q £ Z s and R' C P U R 
such that 

{\p-r\, r£R}={\q-r'\, r' £ R'}. 

As for every pi,p% £ P we have \q— pi\ 2 = \q—p2\ 2 (mod 16), the fact \p— r\\ 2 ^ 
\p — r2\ 2 (mod 16) implies that R' can not be entirely in P. Also, q modulo 
2 differs from (0, ... ,0) in an odd number of coordinates while it differs from 
(1, . . . , 0) in an even number of coordinates or conversely. Therefore \q — p'\ 2 ^ 
\q — r\ 2 (mod 2) if p' £ P and r £ R. This rules out the possibility that R' 
intersects properly both P and R. It follows that R' C R and hence R' = R. 
Note that R U {p} is rigid in the sense that the pairwise distances among the 
points are all different The re constructibility and rigidity of R U {p} gives that 
q = p and hence the claim. 

6 A scheme with shiftable base 

In this section we show that on any finite interval of Z of length n, a simulation 
of an [n/4j-qubit circuit can be done efficiently using global one-qubit gates and 
global two-qubit gates with Hamiltonians of the form A Cj with O < \j\ < 22. 

Here = Z and G = Z acts on fi by translation. Then O = {C 3 \j £ Z}\{C } 
where Cj = {(u, v) £ D x D \ u — v = j}. Let P = 4Z, and for every £ £ 7L let 
Re = {M +1,4^? + 3, A£ + 9}. Just like in the argument preceding Corollary[2j it 
is sufficient to show how to simulate local one-qubit gates at point of P because 
the interaction graph of the global 2-qubit Hamiltonian A Ci connects every pair 
of points in P. 

During simulation of a one-qubit operation at position At, the three qubits 
at positions from Ri are set to one and the qubits at position from Z\ (PL) Re) 
are set to zero. This is a state admitted by the scheme (P,Re). By Proposi- 
tion [T21 (P, Re) is a strictly addressable scheme. In particular, A£ and 4£ + 4 
are 3-addressable. Therefore, by Theorem (TJ any one-qubit gate can be effi- 
ciently approximated using global one-qubit gates and global two-qubit gates 
with Hamiltonians of the form A Cj with O < \j\ < 9. 

After each step of simulation, we have a state admitted by (P, Re) for some 
£ £ Z. (Initially £ = 0.) To perform another simulation step, i.e., to simulate a 
one-qubit gate at position 4£', we have to produce the state admitted by (P, Re>) 
where the state of the qubits at positions from P are left unchanged. Intuitively, 
the configuration Re has to be shifted to Re' . Below we show a procedure doing 
this for £' — £ + 1. For general £' we need \£— £'\ iterations of this basic procedure 
or its inverse. 

In every intermediate step during the shifting procedure the qubits at three 
positions {r±, r2, r%} are set to one while the value of the qubit at a fourth posi- 
tion r 4 is flipped. In each step the system {n, r 2 , r 3 }, satisfies the conditions 
of Lemma [TBI and also |r*j — Tj\ < 22. Therefore the one-qubit operation flipping 
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the value of the qubit at r 4 can be efficiently approximated using global one- 
qubit gates and and global two-qubit gates with Hamiltonians of the form A Cj 
with O < \j\ < 22. Below we give the sequence of consecutive configurations 
{n,r-2,r 3 },r 4 . 

Step 1. {r u r 2 ,r 3 } = {U + 1,4* + 3, At + 9}, r 4 = At + 23 (r 4 : ^ 1). 

(Here {\n - r } \} = {2,6,8} U {22, 20, 14}, n + r 2 + r 3 = 121 + 16 # 121 + 69 = 3r 4 .) 

Step 2. {n, r 2 ,r 3 } = {U + 3,4^ + 9,4^ + 23}, r 4 =4£+l (r 4 : 1 ^ 0). 

(Here {\n - r } \} = {6, 14,20} U {2,8,22}, n + r 3 + r 3 = 12^ + 35 # 12f + 3 = 3r 4 .) 

Step 3. {ri,r 2 ,r 3 } = {4^ + 3,4^ + 9,4^ + 23}, r 4 = At + 5 (r 4 : i-> 1). 

(Here {|r< - r,|} = {6, 14,20} U {2,4, 18}, n + r 2 + r 3 = 121 + 35 # 12f + 15 = 3r 4 .) 

Step 4. {n,r-2,r 3 } = {4£ + 3,4£ + 5,4£ + 23}, r 4 = 4£ + 9 (r 4 :1^0). 

(Here {|r< - r,|} = {2, 18,20} U {6,4, 14}, n + r 2 + r 3 = 121 + 31 # 12f + 27 = 3r 4 .) 

Step 5. {n,r 2) r 3 } = {At + 3, At + 5, At + 23}, r 4 = 4£+ 11 (r 4 : 1). 

(Here {|r< - r 3 -|} = {2, 18,20} U {8,6, 12}, ri + r 2 + r 3 = 12f + 31 # 12£ + 23 = 3r 4 .) 

Step 6. {n,r 2j r 3 } = {4£ + 3,4^ + 5,4^+11}, r 4 = 4£ + 23 (r 4 : 1 h-> 0). 

(Here {|r< - r 3 '|} = {2, 6, 8} U {20, 18, 12}, ri + r 2 + r 3 = 12f + 16 j4 12£ + 69 = 3r 4 .) 

Steps 7-12 are the same as Steps 1-6, respectively, with n + 2 in place of 7*j. 
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